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Abstract 
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206(4)tl^ MeV, /b, = 220(6)ti MeV, = 229(3)1^ MeV and /k = 150(3)^12 
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splitting and D meson masses are also presented. 
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1 Introduction 



The accurate determination of the B and D meson decay constants is of pro- 
found importance in phenomenology. The combination /bV-^b, for both 
and Bg mesons, plays a crucial role in the extraction from experimental data 
of CKM quark mixing and CP violation parameters. The phenomenological 

parameter describes B"-B° mixing, and is expected to be close to unity. 
D meson decay constants are needed for calculations based on factorisation of 
non-leptonic B meson decays to charmed mesons (see [1] for a review and [2] 
for a recent application). 

Numerical simulations of lattice QCD provide a method for computing the 
requisite matrix elements from first principles. A prime concern in such cal- 
culations is the control of discretisation errors, most notably those associated 
with heavy quarks, given that in practice the 6-quark mass is large in lat- 
tice units: amq > 1. One approach, pioneered by Bernard et al. [3] and by 
Gavela et al. [4] and used in the present study, is to work with heavy quarks 
around the charm quark mass and extrapolate to the b mass scale guided by 
continuum heavy-quark effective theory (HQET). Other techniques use some 
form of effective field theory directly to reduce the cut-off effects associated 
with heavy quarks. Examples include lattice HQET [5], non-relativistic QCD 
(NRQCD) [6] and a re-interpretation of the Wilson fermion action as an ef- 
fective theory for heavy quarks, developed by El Khadra et al. [7,8], known as 
the Fermilab formalism. A critical discussion of these various methods, their 
associated systematic errors and a survey of recent results for /b, is given in 
a comprehensive review by Bernard [9] . 

This paper presents decay constants and masses of heavy-light mesons cal- 
culated in the quenched approximation to QCD at two values of the lattice 
coupling, /? = 6.2 and (3 — 6.0. The calculation uses a non-perturbatively 
improved relativistic Sheikholeslami-Wohlert (SW) [10] fermion action and 
current operators so that the leading discretisation errors in lattice matrix 
elements appear at O(a^) rather than 0{a). However, non-perturbative im- 
provement does not necessarily reduce lattice artefacts at a given (3. Given 
results at two /3 values, heavy quark mass-dependent discretisation errors can 
be estimated by combining a continuum extrapolation with the fit in heavy 
quark mass, although a simple continuum extrapolation of the final result is 
not attempted. 

Details of the lattice calculation and the extraction of decay constants and 
masses from Euclidean Green functions are described in Section 2. The ex- 
trapolations to physical quark masses, both heavy and light, and the heavy 
quark symmetry (HQS) relationship between pseudoscalar and vector decay 
constants are discussed in Section 3. Results for the decay constants are pre- 
sented in Section 4 and summarised here. 
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195(6)li MeV 
206(4)1}^ MeV 
220(6)11 MeV 
229(3)lfi MeV 




28(1)11 

8.6(3)1^ 

25(l)li 

8.3(2)li 




1.13(1)1^ 

l-ll(l)l3 

150(3)1^1 MeV 




1.10(2)12 
1.04(1)11 



The first error quoted is statistical, the second systematic. Systematic errors 
are discussed in Section 4, with the main contributions itemized in Table 11. 



2 Details of the calculation 

2.1 Improved action and operators 

In the Wilson formulation of lattice QCD, the fermionic part of the action 
has lattice artifacts of 0{a) (where a is the lattice spacing), while the gauge 
action differs from the continuum Yang-Mills action by terms of 0{a^). To 
leading order in a the Symanzik improvement program for on-shell quantities 
involves adding the SW term to the fermionic Wilson action. 



Full 0{a) improvement of on-shell matrix elements also requires that the cur- 
rents are suitably improved. The improved vector and axial currents are 



(1) 



X 



Al{x) = A^{x) + acAd^P{x) 



(2) 



where 



Vf,{x) =ij{x)jf,ij{x) 

^/.(a;) =^(a;)7^75^(a;) 

P{x) =Tp{x)'~^^1p{x) 

T^y{x)=i){x)ia^yi){x) 
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and dfjt is the symmetric lattice derivative. The generic current renormahsation 
is as follows {J — A,V): 

j''^Zj(l + bjam,)j' (3) 



where Zj is calculated in a mass-independent renormahsation scheme. 
The bare quark mass, aniq, is 



-(- 

2 \K ^crit 



""^9 = o ( ~ ~ Z ) (4) 



where k is the hopping parameter. For non- degenerate currents, an effective 
quark mass is used in the definition of the renormalised current, corresponding 
to 

+ (5) 



In this renormahsation scheme, the improved quark mass, as used in chiral 
extrapolations, is defined as 

mq^mq{l + hmamq) (6) 



2.2 The static limit 

The static quark propagator is calculated using the method of Eichten [5], 
keeping only the leading term in the expansion of the propagator in inverse 
powers of the quark mass: 

5Q(x,i;0,0)= f^(t)e-«*i^ + ^(-t)e™'3*l_jl')5(3)(f)To(t,0) (7) 



where 0) is the product of temporal links from (0, t) to the origin. 

To(i,0)= nt/l(0,i-T) (8) 

r=l 

The prescription for renormalizing and improving the axial static-light current 
is largely analogous to the propagating heavy- light current [11] case. The dif- 
ference is that the temporal static axial current requires covariant derivatives: 

{Af^% = At' + acTi^i^kl,\ (hk + hi) (9) 
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so that 



(Af ^)o = ZT\l + hTam^){Af^)o (10) 

The Ca^* term is not implemented in this calculation, but its value is negative 
which will be significant when the static point is compared to the infinite mass 
extrapolation. 

2.3 Definitions of mesonic decay constants 

The pseudoscalar and vector meson decay constants, fp and /y, are defined 

by 

(0|Aj(0)|P) = ^p,/p, (11) 
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{0\V^mV)^ie,^ (12) 
Jv 

where |P) is a pseudoscalar meson state with momentum p^, while \V) is a 
vector meson, with mass aMy and polarisation vector e^. (V^) denotes 
the renormalised axial (vector) current, here taken to be the renormalised, 
improved local lattice axial (vector) current, defined via equation (3). 



2.4 Simulation details 

In this study, gauge field configurations are generated using a combination of 
the over-relaxed [12,13] and the Cabibbo-Marinari [14] algorithms with pe- 
riodic boundary conditions at two values of the gauge coupling (3 — Q/qq. 
At each (3, heavy quark propagators are computed at four values of the hop- 
ping parameter, corresponding to quarks with masses in the region of the 
charm quark mass. For light quark propagators, three values of k, are used, 
corresponding to masses around that of the strange quark. Table 1 lists the 
input and derived parameters. Errors quoted in this and other tables are sta- 
tistical only unless otherwise specified. The value of the hopping parameter 
corresponding to zero quark IHclSS, ^crit 5 

taken from [15]. The strange and 
normal quark masses are fixed using the pseudoscalar meson masses for the 
pion and kaon where the lattice spacing has been fixed from the pion decay 
constant, /tt. Statistical errors are estimated using the bootstrap [16] with 
1000 re-samphngs. 



2.5 Improvement coefficients 

The improvement program requires values for the current and mass improve- 
ment and renormalisation coefficients defined in equations (2), (3) and (6) in 
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Table 1 

Input and derived parameters. The lattice spacing is set by /jr- 

P = 6.2 /3 = 6.0 

Volume 24^ x 48 16^ x 48 

csw 1-614 1.769 

iVeonfigs 216 305 

(GeV) 2.66l^ 1.91^^ 

Heavy K 0.1200,0.1233,0.1266,0.1299 0.1123,0.1173,0.1223,0.1273 

Lights 0.1346,0.1351,0.1353 0.13344,0.13417,0.13455 

Kcrit 0.135811? 0.135251? 

Kn 0.135771? 0.135191? 

0. 1317911^ 0.13381} 



Section 2.1. One would like to use non-perturbative determinations of the coef- 
ficients in order to remove all 0{a) errors. However, different non-perturbative 
determinations may give mixing coefficients differing by terms of 0{a) and 
normalisation coefficients differing at (9(a^). Using a consistently determined 
set of coefficients (that is, applying the same improvement condition at each 
value of the coupling and a consistent set of improvement conditions for all 
coefficients) should enable a smooth continuum extrapolation. 

Alternatively, perturbation theory may be used [25], although this leaves resid- 
ual discretisation errors of 0{a'^a). Lattice perturbation theory is improved 
by using a boosted coupling g"^ = Qq/uq [26]. The mean link, uq, is related to 
the plaquette expectation value, Uq = (ReTr[/p)/3. 

Non-perturbative determinations of the improvement coefficients are available 
from two groups. The ALPHA collaboration have determined the value of 
Csw [17,18] using chiral symmetry and Ward identities in the Schrodinger 
Functional (SF) formahsm. They have determined ca [18] and Za, Zy and 
by [19] in the same scheme and have a preliminary determination of cy [20,21]. 
Bhattacharya et al. [22-24] have determined all the improvement coefficients 
needed to improve and renormalise quark bilinears, also using Ward identities, 
but on a periodic lattice with standard sources. They also use the ALPHA 
value of Csw to improve the action. 

The two non-perturbative determinations of the improvement coefficients give 
very similar values for the renormalisation coefficients (^'s) and the quark 
mass constant hy (Table 2). However, the mixing coefficients ca and cy differ 
greatly for the two values of the coupling used here. This is shown in Figure 1. 
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Table 2 

Improvement coefficients from various determinations. LANL refers to Bhattacharya 
et al. and BPT is boosted perturbation tlieory. 







p = o.z 






p = o.U 






T A TVTT 


A T U) LJ A 


or 1 


T A ATT 


A T O LJ A 


or 1 


7 A 


v.oLoyz J yO j 


U.OU 1 yO ) J 


U.OIUO 


U.OU / y^ J yo j 


U. 1 c/UU\^c/^ ) 


u.ouoo 


bA 


1.32(3)(4) 




1.24 


1.28(3)(4) 




1.26 


CA 


-0.032(3)(6) 


-0.038(4) 


-0.012 


-0.037(4)(8) 


-0.083(5) 


-0.013 


Zv 


0.7874(4) 


0.7922(4)(9) 


0.7959 


0.770(1) 


0.7809(6) 


0.7820 


bv 


1.42(1)(1) 


1.41(2) 


1.24 


1.52(1) 


1.54(2) 


1.26 


cv 


-0.09(2)(1) 


-0.21(7) 


-0.026 


-0.107(17)(4) 


-0.32(7) 


-0.028 



• ALPHA — PT 

■ Bhattacharya ef al. ♦ BPT 

1* 



6 7 6 7 

P P 

Fig. 1. The improvement coefficients ca (left) and cy (right) as functions of the 
coupling p. Note that the vertical scales in the two plots are different. The values 
of the coefficients are shown for the two values of the coupling used in this work. 
The dashed line through the ALPHA ca points is their interpolating function. 

Changes in ca and cy can have a particularly large effect on the extracted 
values of the decay constants when the pseudoscalar and vector meson masses 
are not small (in lattice units), since the improved current matrix elements 
are given by, 

{0\Al\P) ={0\A,\P) +CAsmh{aMp){0\P\P) 

{0\Vl\V, e) = {0\Vi\V,e) + cy sinh {aMy) (0|T,4| V^, e) (13) 
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Fig. 2. The mixing of the axial current with the pseudoscalar density. P = 6.0, 
kh = 0.1123, kl = 0.13344. The ratio plotted is defined in equation (19). 

when the ground state is isolated. In particular at f3 = 6.0 the heavy-light 
meson mass at the heaviest kappa is somewhat bigger than one such that 
sinh aM ~ 1.3. This is illustrated for the pseudoscalar decay constant in Fig- 
ure 2. The figure shows the ratio (between PA and PP correlators, defined 
in equation (19) and proportional to /p), for several different values of ca at 
/? = 6.0. It can be seen that using the NP value of ca from the ALPHA col- 
laboration decreases R"^ by ~ 20% relative to the case ca = 0. The plot also 
shows the ratio determined using the NP value of ca from Bhattacharya et al. 



At f3 = 6.0 the disagreement between ALPHA and Bhattacharya et al. in 
the mixing coefficients is striking. Bhattacharya et al. try to estimate some 
of the cut-off effects in their determination of ca by looking at the difference 
between two- and three-point derivatives. They also note that since only 
halves between (3 = 6.0 and /3 = 6.2, where disagreement is substantially 
decreased, even higher order cut-off effects are playing a dominant role. Collins 
and Davies [27] have examined the effect of using higher-order derivatives in 
determining ca on the data set used in this work and point out that the 
difference in the value of ca from different orders of derivative disappears once 
the chiral limit is taken. There is clearly a large 0{a) ambiguity in the value 
of Ca- 

Bhattacharya et al. have computed all the coefficients needed for this calcu- 
lation in a consistent fashion and their determinations are used here for that 
reason. The values of the mixing coefficients ca and cy are an order of mag- 
nitude larger in the ALPHA determinations than in BPT at /3 = 6.0. This 
has a large impact on the HQS relation between, and scaling of, the decay 



8 



constants. When the values of the mixing coefficients from Bhattacharya et 
al. are used the HQS relation is satisfied and there is good scaling. 



The ALPHA collaboration compute hy in the SF scheme at very small quark 
masses and in Bhattacharya et al. the heaviest quark mass is amq ~ 0.13. 
One may question whether these coefficients can be used around the charm 
mass: at /? = 6.0, am^ ~ 0.75 and at /? = 6.2, arric ~ 0.5. However, the 
effective normalisation of the vector [28,29] and axial [30] currents has been 
measured at these quark masses and found to be in good agreement with the 
determinations of Zy and hy^. 

A value for hm is required to compute the rescaled quark mass used in the 
chiral extrapolations. Whilst there is a non-perturbative determination [31], 
it is only at one value of the coupling, and so this work follows [15] and uses 
the BPT value with which the non-perturbative value is consistent. 

Expressions for the static improvement and renormalisation coefficients c^*^*, 
6^^* and Z^**^*, have been derived in perturbation theory [11], 



ZT = 1-0 + - 0.137(1)) / + 0{g') (14) 

bT = l~ 0.056{7)gl + O{g',) 
cr = -^xl.01(5)^7o' 

where go is the bare lattice coupling. The one-loop PT and NP renormalisation 
schemes have been matched at the scale /i = rrib [32,33] in the MS scheme. The 
coupling g is then evaluated in the MS scheme [34] at the scale nib. The values 
for foj*^* ^iid cj*^* have been evaluated with the boosted coupling at /5 = 6.0. 
The values are 



Z^*^* = 0.6631}^ (15) 
c5^t = -0.1354 

where the errors on Z^*^* reflect the uncertainty in nib = 4.0 — 4.4 GeV [34] 
and in the MS coupling. 



Indeed Bhattacharya et al. compare their results at /? = 6.2 to [28] 
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Table 3 



Pseudoscalar and vector masses in lattice units at /3 = 6.2 and /3 = 6.0. Fit ranges 



are 12 — 


22 at p - 


= 6.2 and 10-22 (P)^ 


,10-21 (V) at /? = 6.0. 










= 6.2 






/3 = 


6.0 








aMp 


aMv 






aMp 


aMv 




0.1346 


0.84lt} 


1 O 

O.871I2 




0.13344 


I.I45I1 


I.I88I2 


0.1200 


0.1351 


0.823li 


O.856I2 


0.1123 


0.13417 


i.mtl 


1.166l^ 




0.1353 


0.8171^ 


0.84812 




0.13455 


1.1 lOl^ 


1.158l^ 




0.1346 


0.7391} 


0.77511 




0.13344 


1.0061? 


1.056li 


0.1233 


0.1351 


U.Yzll]^ 


0. 759^2 


0.1173 


U. 13417 


0.98ll^ 


1.034l^ 




0.1353 


0.714^1 


0.752l^ 




0.13455 


0.969l^ 


1.026l^ 




0.1346 




U.D r o_2 




U. 13344 


0.8511} 


0.915li 


0.1266 


0.1351 


0.6091} 


0.656l^ 


0.1223 


0.13417 


0.8251? 


0.892l^ 




0.1353 


0.6021^ 


0.650li 




0.13455 


0.81ll^ 


0.883l^ 




0.1346 


0.5051} 


0.563l^ 




0.13344 


0.6751} 


0.759l^ 


0.1299 


0.1351 


0.4841} 


O.546I2 


0.1273 


0.13417 


0.6461? 


0.736l^ 




0.1353 


0.4761} 


O.540I2 




0.13455 


0.6311? 


O.727I3 



2.6 Extraction of masses and decay constants of heavy-light mesons 

The pseudoscalar meson masses are extracted from the asymptotic behaviour 
of the two-point correlation functions, 



cfp{t,P) = Y.^"'''(^pit,mp\o,o)) (16) 

X 

{Zm) ^^^^^^^^^^ _ ^i^~^y-aE,T/2 

2aEp 

where T = 48 and aEp is the energy of the lowest lying meson destroyed by 
the operator Qf, and created by Qp . The superscript S denotes a smeared, 
or spatially extended, interpolating field operator constructed using the gauge 
invariant technique described in [35]. Zp{p) is the overlap of the operator 
with the pseudoscalar state given by Zp{p) — (0|Jlp(0, 0)|P(p)). Fit ranges 
are established by inspection of the time-dependent effective mass. Tabic 3 
shows the best fit masses in each case. A similar procedure is used to extract 
the masses of vector mesons from correlation functions constructed with the 
appropriate vector operators. 
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The decay constants are extracted from the large time behaviour of differ- 
ent two-point correlation functions at zero momentum. The PA correlation 
function is used for the pseudoscalar decay constant: 



^li(^,0) = ^{Ai{t,x)Q'p{Q,0)) (17) 



X 



2clMp 

where A\ is the time component of the improved axial current operator defined 
in equation (2). The superscript L on the correlator denotes a local operator, 
in this case the axial current. = {0\A\\P[0)) is the overlap of the local 
axial operator with the pseudoscalar state, from which the decay constant is 
extracted using equations (3) and (11). 

Extraction of the vector decay constant involves the large time behaviour of 
the VV correlation function: 



C^v{t,0) = EE(^'(^'^)^?(0'0)> (18) 



1^ cosh{aMv{t - r/2))e-«^^^/2 



V 

where Vj is a spatial component of the improved local vector current operator. 
Again, S denotes a smeared or spatially extended interpolating field operator 
and L a local operator. The factor Zy is the overlap of the local vector operator 
with the vector state, Zy = J2r 4(0|^fc 1^(0) ^)): from which the vector decay 
constant can be extracted via equation (12). 

Matrix elements proportional to the decay constants can be extracted from 
ratios of correlation functions. The ratio 

^ % t^^h(«Mp(r/2 - t)) (19) 



is used for the pseudoscalar case and is shown in Figure 2 for various values 
of ca- The vector decay constant is determined using 



t— »00 Zy 

Cyy{t) Z^ 



Results for the decay constants are shown in Table 4. 



(20) 



An alternative method is to fit several correlators simultaneously allowing an 
estimation of the contamination of the ground state signal by excited states. 
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Table 4 

Pseudoscalar and vector decay constants in lattice units. At /3 = 6.2 the fit ranges 



are 14 - 21 (P) and 15 - 23 (V). At /? 
16 - 23 (V). 


= 6.0 the fit ranges 


are 14 - 21 (P) and 




/5 


= 6.2 








= 6.0 










fv 




hl 




fv 


0.1200 


0.1346 
0.1351 
0.1353 


0.08921^ 
0.08471^ 
0.0832+^ 


9.03110 
9.28l}^ 
9.37+}^ 


0.1123 


0.13344 
0.13417 
0.13455 


0.12441}° 
0.11871}^ 
0.1167+}^ 


8.791^^ 
8.941}^ 
8.89+}^ 




0.1346 


0.08731^ 


7.981 8 




0.13344 


0.1217l}g 


7.801 8 


0.1233 


0.1351 
0.1353 


0.0829+1 
0.08141^ 


8.17+^° 

— y 

8.2311J 


0.1173 


0.13417 
0.13455 


0.11621}? 
0.11421}^ 


7.89+}n 

— iu 

7.841}? 




0.1346 


0.08471? 


6.891 6 




0.13344 


0.11701 I 


6.741 I 


0.1266 


0.1351 


0.08051^ 


7.001 ^ 


0.1223 


0.13417 


O.lllSlig 


6.791 I 




0.1353 


0.0790l| 


7.031 8 




0.13455 


0.10971}^ 


6.72li| 


0.1299 


0.1346 
0.1351 


0.08041^ 
0.07651^ 


5.761 i 
5.791 5 


0.1273 


0.13344 
0.13417 


0.11021 I 
0.10531 I 


5.621 I 
5.591 6 




0.1353 


0.07511^ 


5.791 I 




0.13455 


0.1032li? 


5.511 ? 



An eight-parameter fit is made to Cpf,, Cp^ and Cpp, allowing for ground and 
first excited state contributions. Results for the pseudoscalar decay constant 
are consistent with the single ratio fit. The ratio method is used for central 
values in the following but the difference from the multi-exponential fits is 
quoted in Table 14 below as one measure of systematic error. 



2.6.1 The static-light axial current 

The static-light axial current can be extracted from axial-axial correlation 
functions. The local-local and local-smeared correlation functions were used 
in this work, with large-time behaviour given by 



C'Ai{t,0) = ^{A^,{t,x)A',\0,0)) (Zi)^e-«^^* (21) 

X 
X 
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LU 

< 
CO 



1.4 - 

1.2 - 

1 - 

0.8 - 

0.6 - 

0.4 - 

0.2 - 









• LL 

o LS 
-■ LLfit 
- LSfit 
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14 
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Fig. 3. The effective mass plot for the LL and LS static-light axial correlation 
function, kl = 0.13344. The lines on the plot show the value of the fit, and the 
fit ranges. 

where AE is the unphysical difference between the mass of the meson and the 
mass of the bare heavy quark. The static-light amplitude, Z^, is 



Stat, 



(22) 



The smearing function is again the one described in [35]. 

The static-light correlation functions were generated at only one value of the 
coupling, (3 = 6.0, and without the covariant derivative operators necessary to 
improve the current. The assumption that the matrix element of the improve- 
ment term is of the same order of magnitude as the primary term leads to 
(9(10%) error in the static point associated with the absence of improvement. 

The signal of the static-axial current is very quickly overwhelmed by statistical 
noise, making the fit difficult. Simultaneously fitting to both the local-local 
and local-smeared correlator gives a better estimate of the local current. This 
is shown in Figure 3. The renormalised values of Zi (without the current 
improvement term) are shown in Table 5. 



2. 7 Extraction of masses and decay constants of light-light mesons 

The masses and decay constants of the light-light pseudoscalar mesons are 
extracted from a simultaneous fit to three correlation functions: Cpf,, Cpp 
and Cp\. The smearing function used for the light-lights is the 'fuzzing' of 
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Table 5 

Values of aAE and Zl from simultaneous fits to the LL and LS correlation functions. 
The fit ranges are 8-12 for the LL and 4-12 for the LS. 





Zl 


aAE 


0.13344 


0.30(2) 


0.80(2) 


0.13417 


0.28(2) 


0.78(3) 


0.13155 


0.2G(3) 


0.7(3(3) 



Tabic 6 

The light-light pseudoscalar decay constants. The fit ranges are 8 — 22 at both values 
of the coupling, except for the lightest mass combination at /3 = 6.0 where the fit 
range was 7—23. 





/3 = 6.2 






/3 = 6.0 








afp 






afp 


0.13460 
0.13510 


0.13460 
0.13460 


0.06401^ 
0.0610i^ 


0.13344 
0.13417 


0.13344 
0.13344 


0.0841^^° 


0.13530 
0.13510 


0.13460 
0.13510 


0.05991^ 
0.05871^ 


0.13455 
0.13417 


0.13344 
0.13417 


0.08211^° 
0.0798li^ 


0.13530 


0.13510 


0.05751^ 


0.13455 


0.13417 


0.0775li| 


0.13530 


0.13530 


0.05591^ 


0.13455 


0.13455 


0.0730;i^ 



reference [36]. The results for the decay constants are shown in Table 6. The 
results for the pseudoscalar masses are compatible with those in [15] where a 
similar fit was used. 



3 Extrapolation and interpolation in the quark masses 

Extrapolation or interpolation in quark masses must be performed to ex- 
tract physical masses and decay constants. For heavy quarks, the presence 
of O(a^mQ) lattice artefacts when using the SW action with the NP improved 
renormalisation scheme imposes a constraint amg < 1 on the quark masses 
that can be studied. This limits hadron masses to ~ 2 GeV or slightly greater 
at f3 = 6.0, see Table 10. Input light quark masses are kept above ms/2 to 
avoid critical slowing down of quark propagator calculations and the possible 
appearance of finite volume effects. Interpolations to rric and rris are needed, 
together with extrapolations for mb and the light quarks. 



14 



Table 7 

The inverse lattice spacing, (GeV) 

P = 6.2 13 = 6.0 

ro 2mtl 2.12ti 

rup 2.54l^ 1.89t| 



3.1 The light-light sector 



The dependence of the hght pseudoscalar decay constant on the hght pseu- 
doscalar meson mass is described by chiral perturbation theory, giving the 
following ansatz: 

afp = Co + Ci{aMpf + C2iaMpy (23) 



The value of afp for which 

aff 




(24) 

expt 



can then be determined and compared to the experimental value of /^r, taken 
from [34], to fix the lattice spacing. At both values of the coupling, quadratic 
and linear fits in {aMpY give essentially the same answer for the lattice spac- 
ing. Thus the data satisfy lowest-order chiral perturbation theory. Because of 
quenching and other systematic effects, determinations of the lattice spacing 
from different quantities disagree. In this work the lattice spacing is set by /tt. 
Lattice spacings fixed by the Sommer scale, Tq, [37,38], and rUp as determined 
by [15] are used to estimate systematic error from this source. The values are 
shown in Table 7. 

The light, or 'normal', quark mass mn, defined by = {m^ + md)/2, and 
strange quark mass rus are determined using the lowest-order chiral perturba- 
tion theory relation for the mass of a pseudoscalar meson with quark content 
qi and g2, 

{am-pY = B{amq^ + arhq^) (25) 



where the rescaled quark mass mq is defined in equation (6). The value of /tent 
is taken from a previous UKQCD calculation [15] and is listed in Table 1. The 
value of the hopping parameter corresponding to the normal quark mass is set 
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0.065 



0.06 
0.055 



0.05 

0.01 0.02 0.03 

Fig. 4. Dependence of the light pseudoscalar decay constant fp on quark masses. 
The hnes give fp as a function of one of the quark masses, arhq^, for fixed values of 
the other quark mass, afhq^. Filled symbols are the lattice data points. The lowest 
line shows the extrapolation of arrig^ to the normal quark mass, with the open 
triangle giving /k from the position of the kaon on this line. 

by the charged pion according to 

a^(m^)expt = 2Bam^ (26) 

and that for the strange quark mass by 

a\ml)j = B{am^ + am,) (27) 

where 

(m^), = ^(m^±+m^o) (28) 

with the lattice spacing set by f-„. These hopping parameters are also listed 
in Table 1. 

The ansatz for the dependence of the light pseudoscalar decay constant on 
quark masses is: 

afp = fo + /i(«^<?i + arhg^) (29) 
Figure 4 shows the results of fitting to this ansatz and extrapolating to /k- 



K, =0.1 3460 




---- K,=0. 13510 


- 


K^=U.1oOoU 




--- ^^=\ 




A K2=Kj, (Kaon) 




^^^^ + 
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0.015 0.02 0.025 0.03 0.035 



0.09 L 

0.085 

0.08 

0.075 
0.84 

^ 0.82 
CO 

0.8 

0.005 0.01 0.015 0.02 0.025 0.03 0.035 

Fig. 5. The chiral extrapolation of the pseudoscalar decay constant (top) and pseu- 
doscalar mass (bottom) against rescaled Hght quark mass. P = 6.2 and kh = 0.1200. 
The vertical lines show the strange and normal quark masses. 

3.2 Chiral extrapolations of heavy-light masses and decay constants 

A linear dependence of heavy-light masses and decay constants on the light 
quark mass is assumed: 

aKi = ai + PiOmq (30) 

where Ki is fp, /y/a, mp or my. Some example extrapolations are shown 
in Figure 5. The results for the extrapolated quantities are shown in Ta- 
bles 8 and 9. 



3.3 Heavy quark extrapolations 

Heavy Quark Symmetry (HQS) implies asymptotic scaling laws [39] for the 
decay constants in the infinite heavy quark mass limit. Away from this limit, 
heavy quark effective theory ideas motivate the following ansatze for the de- 
pendence on the heavy meson masses: 





^p{Mp 
$y(My; 



= e(MB, Mp)fp^/Mp = 7p h + + 



Mf 



Ml, 



V 



Ml 



V 



(31) 
(32) 
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Table 8 

Masses and decay constants at physical light quark masses. /3 = 6.2. 



kh k,l aMp afp aMy fv 



0.1200 




0.800(2) 


0.079(1) 


0.832(4) 


9.6(2) 


Kg 


0.828(2) 


0.0876(8) 


0.860(2) 


9.1(1) 


0.1233 




0.696(2) 


0.077(1) 


0.736(3) 


8.4(1) 


Kg 


0.726(1) 


0.0857(8) 


0.764(2) 


8.05(9) 


0.1266 




0.583(2) 


0.075(1) 


0.634(3) 


7.1(1) 


Kg 


0.615(1) 


0.0832(7) 


0.661(2) 


6.93(7) 


0.1299 


l^n 


0.455(2) 


0.071(1) 


0.523(3) 


5.8(1) 


Kg 


0.491(1) 


0.0789(7) 


0.551(1) 


5.77(6) 



Table 9 

Masses and decay constants at physical light quark masses. /3 = 6.0. 







aMp 


afp 


aMy 


fv 


0.1123 


l^n 


1.087(4) 


0.112(2) 


1.138(5) 


9.0(2) 


Kg 


1.125(2) 


0.122(1) 


1.171(3) 


8.8(1) 


0.1173 


l^n 


0.945(3) 


0.109(2) 


1.005(5) 


7.9(2) 


Kg 


0.985(2) 


0.119(1) 


1.039(3) 


7.8(1) 


0.1223 




0.786(3) 


0.105(1) 


0.862(5) 


6.8(1) 


Kg 


0.829(2) 


0.115(1) 


0.897(3) 


6.75(8) 


0.1273 




0.603(2) 


0.099(1) 


0.704(5) 


5.5(1) 




0.651(2) 


0.1080(9) 


0.710(3) 


5.60(7) 



where © denotes logarithmic corrections given at leading order by [40], 



Here, f3o is the one-loop QCD beta function coefficient, equal to 11 in the 
quenched approximation, and = 295 MeV [41]. 

HQS also relates the pseudoscalar and vector decay constants as follows [40] ; 

Um^f^=(l + l^ + 0(l/M)) (34) 
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where M = (Afp + 3Mv)/4: is the spin-averaged heavy meson mass. The one- 
loop factor in equations (31) and (32) cancels in the ratio in equation (34). 
Higher-order QCD corrections produce the term proportional to ag. U (M) is 
defined to ehminate the radiative corrections in U{M), 

U{M)^U{M)/^l + l^^ (35) 



Calculated values of U{M) are fitted to the following parameterisation: 

U(M) = a;o + ^ + ^ (36) 
M M 



HQS implies that ujq = 1. However, uq can also be left as a free parameter to 
test the applicability of HQS. Likewise, HQS can be applied to set 7p = 7v' 
in fits using equations (31) and (32). However, higher order QCD corrections 
would modify this in a similar way to U (M), that is 

7y V 3 47r y ^ ' 



The systematic error arising from extrapolating the decay constants to the h 
mass is studied in Section 4. The error from truncating the expansion in inverse 
powers of the heavy mass is considered and the issue of the propagation of 
discretisation effects under this extrapolation is addressed. 



4 Decay constants 

The main results for the decay constants are listed in Table 11 and summarised 
in the introduction. Central values are obtained (3 — 6.2, setting the scale 
with /tt. These results are discussed in more detail here, starting with /k. 

For /k the largest source of uncertainty is the choice of quantity used to set 
the scale. The value of /k increases by 8% when the scale is set by ro and 
decreases by 4% when it is set by rUp. When is used to fix the strange 
quark mass rather than mx, the value of /k decreases by 2.5%. With only two 
values of the lattice spacing a continuum extrapolation is not attempted. The 
decay constant is smaller at /? = 6.0 by 2.7%, which is taken as an estimate 
of its discretisation error. These estimates of systematic errors are combined 
in quadrature. 

The ALPHA collaboration [42] have computed the decay constant, /k, with 
the same action as used here, at several values of the coupling. They perform 
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0.44 



0.42 




• UKQCD 
♦ ALPHA 
■ expt 

o ALPHA CL 



0.38 



0.36 







0.01 




0.03 



0.04 



Fig. 6. The dependence of fx on the lattice spacing. 



a continuum extrapolation of tq/k against a^/rg. To examine the scaling be- 
haviour of the decay constant in this work, a comparison with the ALPHA 
results is shown in Figure 6. The line shows the linear extrapolation to the 
continuum limit (CL) that was performed by ALPHA, excluding the point at 
the coarsest lattice spacing, f3 = 6.0. There are two main differences in the 
calculations. First, ALPHA use degenerate light quark masses, whereas this 
work specifically takes into account the non-degeneracy of the quarks. How- 
ever, ALPHA have checked on the same dataset analysised in this work that 
using degenerate light quarks has negligible effect. A more obvious difference 
is in the values of the improvement coefficients, ca and Ba- As noted earlier, 
improvement coefficients may differ by 0{a) terms, depending on the improve- 
ment condition used. A particular choice of conditions used to determine the 
coefficients may result, for a particular quantity, in smaller discretisation ef- 
fects at finite lattice spacing. This seems to be the case for the decay constants 
when the improvement coefficients of Bhattacharya et al. are used rather than 
those determined by ALPHA, especially at /3 = 6.0. However, as the con- 
tinuum limit is approached the values of /k from the two studies converge. 

The major source of uncertainty in the heavy-light decay constants comes 
from the determination of the lattice spacing. Additional systematic errors 
arise from discretisation effects and from the heavy quark extrapolations. Since 
these latter errors are related to the size of the heavy quark masses, the heavy 
masses used in this work are listed in Table 10. All these errors will now be 
discussed. 

The ambiguity in the decay constants from the lattice spacing is quite large: 
there is an overall 16% variation in the value of /b when the scale is set by rrip. 
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Table 10 

Heavy quark and meson masses used in this work. The table shows the bare quark 
mass in lattice units, the renormalisation group invariant quark mass defined by 
= Zjn{l + bmamQ)mQ[43] in GeV and the heavy-light («;„) pseudoscalar meson 
mass in GeV. The scale is set by f^^. 





(3 = 


6.2 






(3 = 


6.0 




K 


arriQ 




Mp 


K 


aniQ 




Mp 


0.1200 


0.485 


1.59 


2.128(5) 


0.1123 


0.756 


1.26 


2.076(8) 


0.1233 


0.374 


1.36 


1.851(5) 


0.1173 


0.566 


1.18 


1.805(6) 


0.1266 


0.268 


1.06 


1.551(5) 


0.1223 


0.392 


0.97 


1.501(6) 


0.1299 


0.168 


0.72 


1.210(5) 


0.1273 


0.231 


0.65 


1.152(4) 



Table 11 

The decay constants at both /? values with the scale set by /tdTq and rup and a 
quadratic (Q) or linear (L) fit for the extrapolation in inverse heavy meson mass. 
The quadratic fit uses all four heavy masses, while the linear fit uses the heaviest 
three. The central values are shown in bold face. 







P = 


6.2 






P = 


6.0 






/p(MeV) 


fv 


/p(MeV) 


fv 






Q 


L 


Q 


L 


Q 


L 


Q 


L 


B 


195(6) 


178(5) 


28(1) 


29(1) 


200(8) 


182(7) 


26(2) 


28(2) 


D 

B, 


206(4) 


207(4) 


8.6(3) 


8.6(3) 


210(4) 


211(5) 


8.3(3) 


8.3(3) 


220(6) 


201(5) 


25.2(8) 


26.5(8) 


222(7) 


202(7) 


25(1) 


26(1) 


Ds 


229(3) 


230(4) 


8.3(2) 


8.3(2) 


230(4) 


230(4) 


8.3(2) 


8.3(2) 


B 


217(5) 


200(4) 


24.8(8) 


25.9(7) 


226(6) 


207(5) 


18.9(8) 


20.4(8) 


D 


222(3) 


223(3) 


7.8(2) 


7.8(2) 


229(3) 


230(3) 


6.7(2) 


6.8(2) 


Bs 


240(4) 


222(3) 


22.9(5) 


24.0(4) 


247(4) 


226(3) 


21.9(7) 


23.1(5) 


Ds 


243(2) 


244(2) 


7.7(1) 


7.7(1) 


246(2) 


247(2) 


7.5(1) 


7.5(1) 


B 


185(4) 


168(3) 


29.5(9) 


31.1(9) 


197(6) 


179(4) 


26(1) 


28(1) 


D 


199(3) 


199(3) 


9.1(2) 


9.1(2) 


208(3) 


209(3) 


8.4(2) 


8.4(2) 


rrip 

B, 


209(3) 


191(2) 


26.5(5) 


28.0(4) 


220(3) 


200(2) 


24.9(7) 


26.5(5) 


Ds 


221(2) 


222(2) 


8.7(1) 


8.7(1) 


228(2) 


228(2) 


8.4(1) 


8.4(1) 



ro or (see Table 11). The difference between the decay constants determined 
at /5 = 6.2 and (3 = 6.0 is smallest when the scale is set by as shown in 
Figure 7. This is not surprising since some systematic uncertainties may cancel 
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Fig. 7. The dependence of the pseudoscalar decay constants on the lattice spacing. 
The scale is set by 

in ratios of decay constants. Therefore the central values are quoted using f^, 
to set the scale. 

Discretisation errors are considered next. Although this calculation is 0{a) 
improved, C(a^) effects may still be important at any fixed lattice spacing. In 
particular, (9(a^m^) effects could well be significant for heavy quarks. With 
results at only two values of the lattice spacing, a full continuum extrapolation 
cannot be attempted. The heavy-light decay constant at the finer lattice spac- 
ing is taken as the central value, with the result at the coarser spacing used 
as one indicator of discretisation errors. A fuller investigation of discretisation 
effects for the heavy quarks follows. 

The general characteristic size of discretisation effects in this simulation is 
investigated by examining the free particle dispersion relation, which is altered 
in discrete spacetime. The lattice dispersion relation can be written 



where Mi is the energy at zero momentum and M2 is the kinetic mass, de- 
fined by Ma"^ = d'^E/dpf\^^Q. The dispersion relation is investigated by fit- 
ting hadronic correlators computed at five different momentum values (p ^ = 
0, 1, 2, 3, 4 in lattice units) using linear and quadratic fits in p"^. The fit for the 
heaviest quark combination at /3 = 6.0 is shown on the left of Figure 8. For 
lighter mass combinations at /3 = 6.0 the signal for highest momentum is very 
poor, and so in general the fits exclude the = 4 channel. The right hand 
side of Figure 8 shows the M1/M2 estimates for each heavy quark. Also shown 



(38) 
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Fig. 8. Left: the pseudoscalar dispersion relation for kh = 0.1123, hl = 0.13344 in 
lattice units at /? = 6.0. Right: the ratio M1/M2 as a function of Mi. 



on the right in Figure 8 is the ratio of M1/M2 when M2 has been determined 
by the shift in quark mass from uii to m2 [44]. 

^ = Ml + (m2 - mi) (39) 



and the quark masses mi and m2 have been calculated from the tree-level 
relation described in [7,45,46]. 

ami = + amq) (40) 
' +^ (41) 



am2 amq {2 + amq) 1 + amq 



Here M denotes hadron mass and m denotes quark mass. There is good agree- 
ment between the tree-level perturbative description of M2 and the NP de- 
termination. At the quark masses used in this work the ratio M1/M2 remains 
close to unity, so that deviations from the continuum dispersion relation are 
small. Finally, it is worth noting that M2 is determined at non-zero momen- 
tum, which is statistically noiser, and so has statistical errors an order of 
magnitude larger than for Mi. 

The axial current in this work is normalised according to equations (2,3). 
As already mentioned the improvement term is proportional to sinh(aMp) 
and the mass-dependent normalisation is proportional to amq. Bernard [9] 
proposed the following alternative normalisation, labelled CB, 

(^0 )CB = Za 



1 + bAdm + 2acA 



do{0\P\P) 
rh{0\Ao\P) 



An 



(42) 



where m is the average of the bare quark masses in the heavy-light state. The 
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Fig. 9. The left plot shows the extrapolations of the pseudoscalar decay constant in 
heavy meson mass at /? = 6.0 with different normalisations and masses. The right 
plot shows a comparison between the extrapolations at both values of the coupling. 

CB normalisation differs from the NP improvement scheme normalisation only 
at 0(a^). Taking into account an additional normalisation of y/AnqHg, the CB 
norm has a finite static limit since KQarriQ ^ 1/2 in as kq —>■ 0. The CB norm 
looks very like the normalisation of Kronfeld, Lepage and Mackenzie (KLM 
norm) [7,26,45,46], 

— ^ = ip^l + nam (43) 

where /z would be the mean link in the tadpole improved tree-level Fermi- 
lab formalism. In this case /i is a mass-dependent factor that contains the 
information about improving the current. 

Figure 9 shows the mass dependence of the function $, defined in equation 
(31), for a variety of methodologies. The NP norm is shown with both a 
quadratic fit (Q) to all four masses and a linear (L) fit to the heaviest three. 
The CB norm data has the extrapolation using the M2 masses rather than 
Ml. It is clear from the figure that the difference between the "NP Q vs Mi" 
and "CB Q vs M2" fits is smaller than the difference between the quadratic 
and linear fits to NP norm. It should be noted that the effect of changing the 
normalisation to CB differs from that described in [9] which used the ALPHA 
determination of ca and the preliminary NP determination by Bhattacharya 
et al. [22] of 6^- 

In the 0{a) improved theory used here, 0{a^mQ) lattice artefacts could still 
be large: the worst case simulated has a'^rriQ ~ 0.5 (see Table 10) and one 
may worry that mass-dependent discretisation effects, which may not be too 
large at the charm scale, are enhanced by extrapolation to the bottom scale. 
Although a Taylor expansion of {arnqY in l/mq around the charm quark 
mass needs large terms to reproduce {arribY, the issue is not directly the fate 
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of 0{a'^mQ) under extrapolation, but rather the effect of such terms at the 
charm scale on the extrapolation. 

The left hand plot in Figure 9 includes the static point. The smaller set of 
error bars shows the error from renormalisation, the larger, the statistical 
uncertainty. The static point is higher than the extrapolations; however, the 
missing static improvement term would lower the static point, since c^^^ is 
negative. For this reason an interpolation to the h quark mass using the static 
and c quark data is not implemented; rather the static point is used as a 
check of the extrapolation. The result supports the view that, the increase in 
discretisation error after the heavy quark extrapolation is moderate. 

The right-hand side of Figure 9 shows the extrapolation of $p for both values 
of /3 in physical units. Although the difference between the curves does in- 
deed grow during the extrapolation, it is smaller than the statistical errors at 
the bottom scale. The qualitative agreement between the two curves suggests 
that discretisation errors are not large. To study this further, equation (31) is 
modified as follows, 

$(M, a)=7fl + 4 + ^+ <^^f + C{aMf\ (44) 



and data from both lattice spacings fitted simultaneously. This uses the two 
values of a to increase the number of aM values available, allowing the im- 
portant class of heavy-mass dependent discretisation errors to be studied. The 
sum of the first three terms in equation (44) has these errors subtracted out 
and is hereafter referred to as the 'quasi-continuum' result. Figure 10 shows 
the quasi-continuum curve as a function of 1/M, together with $(M, a) at 
each j3 value separately. The curves A and B which blow up as 1/M 
are the fits to equation (44) containing lattice artefacts. However, the quasi- 
continuum result, where these artefacts are subtracted, does not differ greatly 
from the extrapolations at the individual lattice spacings using equation (31). 
The quasi-continuum $ gives Jb = 186(10) MeV, which differs by 5% at 
(3 = 6.2 (the number appearing in Table 14) and 7% aX j3 = 6.0 from the 
extrapolations using (31). Excluding the C term, or using amq instead of aM 
makes little difference. Using [aMYK/M instead of (aM)^ also makes no sig- 
nificant difference. Varying all five fit parameters in equation (44) over the 
region where the chi-squared per degree of freedom increases by up to 1 from 
its minimum value gives a variation of ±16 MeV for the quasi-continuum /s, 
so the fit is stable. 

The systematic error associated with truncating the HQET series is now dis- 
cussed. First, linear and quadratic fits to the data are compared. The linear 
fit uses the heaviest three quark masses and drops the 7] term from equa- 
tions (31) and (32). This is compared to a quadratic fit using all four masses. 
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Fig. 10. Effects of lattice artefacts on the heavy quark extrapolation. The solid 
line shows the quasi-continuum extrapolation described in the text. The dashed 
and dotted lines show equation (31), when extrapolating at each value of (3 
separately. Line A(B) shows the function $(M, a), equation (44), at P = 6.2(6.0). 



The resulting variations, which are shown in Table 11, amount to a relative 
error of 9% for /b, indicating that a quadratic term is necessary. A cubic fit is 
performed to check that a quadratic term is sufficient. There are not enough 
data points for a cubic fit at each lattice spacing, so the fit combines data from 
both lattices as before, adding a cubic term ^/M^ to $(M, a). This results in 
a value of /b which is 5% greater than the quasi-continuum determination. 
This is also included as a systematic error in Table 14. 

The HQS relation between the pseudoscalar and vector decay constants can 
also be investigated. The quantity U (M) defined in equation (35) should be 
equal to unity in the static limit. The extrapolation of U{M) is shown in 
Figure 11 and displayed in Table 12. "Q" denotes a quadratic fit to all heavy 
quarks whereas "L" denotes a linear fit to the heaviest three. The extrapolation 
for P = 6.0 displays the expected static limit, albeit with large errors. At (3 = 
6.2, U{M) from the quadratic fit deviates from unity in the static limit. On the 
finer lattice, discretisation errors are smaller and so one might expect better 
agreement in the static limit. However, U{M) depends on the ratio of the 
axial to vector currents and is particularly sensitive to the mixing coefficients, 
Ca and Cy, which are the most poorly known of the improvement coefficients. 
Different determinations of cy vary by an order of magnitude, even at jS = 6.2, 
where all the other coefficients are in much better agreement. Indeed, varying 
cy by the quoted errors of 25% in Bhattacharya et al. varies the value of 
[/(Moo) by around 6%. Taking into account the uncertainty in the value of 
cy, the value of [/(Moo) is consistent with unity. 



26 



1.3 

1.2 

1.1 

^ 1 
1^ 

0.9 
0.8 
0.7 
0.6 

0.1 0.2 0.3_ 0.4 0,5 0.6 0.7 0.8 

1/M (GeV) 

Fig. 11. The quantity U{M) as a function of the inverse spin averaged mass M. 
Table 12 

U{M) as a function of heavy meson mass. 







Md 


Mb 


Moo 


6.2 


Q 


0.83(1) 


1.01(2) 


1.13(4) 


L 


0.83(1) 


0.98(2) 


1.05(5) 


6.0 


Q 


0.82(2) 


0.96(4) 


1.04(7) 


L 


0.82(2) 


0.94(3) 


1.01(4) 



Table 13 

Flavour-breaking ratios. 



(3 




/Bs 

/b 


/d 


/bJ 


/dj 


6.2 


(i) 


1.13(1) 


1.11(1) 


1.10(2) 


1.04(1) 


(ii) 


1.11(2) 


1.11(1) 


1.10(2) 


1.05(1) 


6.0 


(i) 


1.11(2) 


1.09(1) 


1.05(3) 


1.00(2) 


(ii) 


1.10(2) 


1.09(1) 


1.05(3) 


1.01(2) 



The SU (3) flavour breaking ratios are determined by two methods: from the 
ratio of extrapolated values of the decay constants, labelled (i) in Table 13, 
or by extrapolating the ratios constructed at each heavy mass, labelled (ii) in 
the Table. In the ratio of fpjfp^ most of the heavy quark mass dependence 
seems to cancel. The fact that the two methods agree, within statistical errors, 
suggests that the systematic error in these ratios is small. 
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Table 14 

Percentage systematic uncertainties. Systematic differences are obtained by varying 
the procedure used to calculate the decay constants. The central values are fixed by 
using to set the scale at /3 = 6.2, performing a quadratic heavy quark extrapo- 
lation, taking the central values of the mixing coefficients c^y, fitting correlation 
functions to a single exponential, and using mp^ to set the strange quark mass. 



Pseudoscalar 


/b 


/d 


/Bs 




/b 


J? 

/d 


scale set by ro 


11 


8 


9 


6 


-2 


-2 


scale set by rup 


-5 


-3 


-5 


-3 


1 


1 


linear vs quadratic 


-9 




-9 








quasi-continuum 


-5 


-2 


-5 


-2 






quasi-continuum plus 


5 




5 








(3 = 6.0 


3 


2 


1 


1 


-1 


-2 


multi-exp 


-2 


-3 


-5 


-3 


-3 


1 


strange quark mass from 






-2 


-2 


-2 


-2 


coeff. CA 


1 


1 


1 


1 


1 




Vector 


/b* 


/d* 


/b* 


/d* 


Jb* 
/b* 


/D* 

/d- 


scale set by ro 


-11 


-10 


-9 


-8 


2 


2 


scale set by rUp 


6 


5 


5 


5 


1 


1 


/3 = 6.0 


-7 


-4 


-2 


-1 


-5 


-4 


linear vs quadratic 


5 




5 








quasi-continuum 


6 


3 


6 


3 






quasi-continuum plus 1/M^ 


-5 




-5 








strange quark mass 






2 


2 


-3 


-2 


coeff cy 


3 


1 


3 


1 


1 


1 



The systematic variations of the decay constants are shown in Table 14, in or- 
der of importance. The total systematic uncertainty is obtained by combining 
these in quadrature. 

4-1 Comparison with other determinations 

A comparison with some other recent quenched results for /b is shown in 
Table 15. For a comprehensive review, the reader is referred to the article by 
Bernard [9]. The results which use the Fermilab formalism [7] are continuum 
limits. The APE result uses the same action at /3 = 6.2 and the UKQCD 
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Table 15 

Comparison with some other recent determinations of /b- "FNAL" denotes use 
of the Fermilab formahsm [7] and "charm" indicates that the heavy quarks have 
masses around charm and are then extrapolated to the bottom scale. NB: the World 
Average is not the average of the numbers displayed below. 



Heavy Quark /b (MeV) 



El-Khadra et al. [47] 


FNAL 




MILC [48,49] 


FNAL 




JLQCD [50] 


FNAL 


173(4)(9)(9) 


UKQCD Tad [51] 


charm 


177(17)^22 


APE [52] 


charm 


179(18)1^^ 


CP-PACS [53] 


FNAL 


188(3) (9) 


This Work 


cllMlMll 


195(6);^^ 



World Average [9] 175(20) 



Tad result uses the same gauge configurations as this work but with a tadpole 
improved value of csw- 

The value of /b obtained in this work is the highest, although it is compatible 
with the world average. One reason for this might be the heavy quark method 
employed. However, other results using the same heavy quark extrapolation 
have a lower value, while the most recent result using the Fermilab method 
from the CP-PACS collaboration [53] has a value in agreement within sta- 
tistical errors. Different heavy quark methods have different associated sys- 
tematic uncertainties. A detailed analysis of those present in this work has 
been discussed above, including an estimate of the effects of the heavy quark 
extrapolation. These uncertainties arc also explored in [9], which additionally 
addresses some of the systematic issues affecting the Fermilab formalism. The 
value obtained in this work is consistent within these systematic uncertainties 
with other quenched results. 

There have been several recent calculations of /b with two flavours [Np = 2) 
of dynamical fermions using NRQCD or the Fermilab formalism. Again, they 
are reviewed in [9] . The effect of unquenching is found to increase the value of 
/b by 10 - 15%. 

The only heavy-light decay constant to have been measured experimentally 

so far is /d, [54] , 

/d, = 260 ± 19 ± 32 MeV (45) 
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Table 16 

The hopping parameter corresponding to the charm quark mass. The label FIT 
refers to the ratio e/A in equation (46) having a fitted value, whilst BPT refers to 
the ratio having the BPT value. 





(3 = 6.2 


P = 6.0 




FIT BPT 


FIT BPT 



a 


-Hro) 


0. 124981^ 


0.12513t^ 


0.119521^ 


0.12056 


+8 
-6 


a' 




0.1222X1^ 


0.122271^ 


0.1160li 


0.1165^ 


h2 
'2 


a' 




0.12321^ 


0.12331^ 


0.11641^° 


0.11711 


3 
10 



With such large uncertainties, the experimental value is consistent with both 
the unquenched and quenched values of /d^. 



5 Spectroscopic quantities 

To determine the value of the hopping parameter, k^, corresponding to the 
charm quark mass requires an interpolation in heavy quark mass. The im- 
proved quark mass definition in equation (6) also applies to heavy quarks, so 
if the meson mass depends linearly on the improved heavy quark mass, the 
corresponding hare quark mass dependence is, 

aM-ii{mQ) = p + XaniQ + ea^mq (46) 



where e/A = bm- In the NP improved formulation all lattice artefacts of 0{a) 
have been removed. However, 0{a'^mQ) effects for the heaviest quarks could be 
significant. Any contributions at this order would affect the ratio e/A such that 
it was no longer equal to 6^- A fit to equation (46) was tried with e/A fixed to 
bm from boosted perturbation theory, and with e/A allowed to vary freely[3. 
The heavy quark dependence can be used to fix by choosing a particular 
state (or splitting) to have its physical value. Choosing the pseudoscalar mass 
to fix Kc, the spectrum of heavy-light mesons can then be predicted. 

For (3 = 6.2, the results of using equation (46) are shown in Figure 12. Whilst 
the value of e/A = —0.505(4) (labelled FIT in the figure) differs somewhat 
from the value of bm = —0.652 from BPT, it makes little difference to the 
value of the Kc, of order 0.1%. The choice of quantity to set the lattice spacing 
clearly has a rather large effect. Values for are shown in Table 16, using the 
free fit to set the value of the hopping parameter. 

^ This procedure is entirely equivalent to that of Becirevic et al. [52]. 
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For (5 = 6.0, the fitted value of e/A = —0.384(3) differs significantly from 
the value of hm = —0.662 from BPT. Consequently the fit using the BPT 
has a very large x^; as the data cannot accomodate a model with such large 
curvature. However, the difference in is still small (of order 1%). For (3 = 6.2 
the heaviest quark has a value of arriQ = 0.485 whereas for (3 = 6.0 the heaviest 
quark has amg = 0.775. Discretisation errors of (9(a^m^) could be responsible 
for modifying the value of the ratio e/A. The value of Kc is rather insensitive 
to e/A but because the value of amq is so large, the improved quark mass 
changes dramatically. Using the free fit as the preferred method, the value of 
Kc is shown in Table 16. 

The spectrum of heavy-light charm states is shown in Figure 13 and tabulated 
in Table 17. With only two values of /3 a continuum extrapolation is not at- 
tempted, but the difference between the two couplings can be used to estimate 
systematic errors. It is clear from the figure that a large systematic error arises 
from the definition of the lattice spacing. Depending on which quantity is cho- 
sen, different mass splittings appear to be closer to their continuum values. 
The central values for the spectrum are produced by the following procedures: 
single exponential fits, using tq to set the scale, with e/A a free parameter 
in equation (46). The values of the masses at /3 = 6.2 are taken to be the 
central values, with the values at /5 = 6.0 used to estimate systematic error. 
The remaining estimates of systematic error come from using multiple expo- 
nential fits, using rUp to set the scale and, in the case of states with a strange 
quark, using mx* instead of m|- to set the strange quark mass. All these sys- 
tematic errors are combined in quadrature. The results, with experimental 
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Fig. 13. The scaling of the D meson spectrum. The filled symbols have the lattice 
spacing set by ro, and the open symbols by nip. 

comparisons, are: 



Lattice, this work 



Experiment [34] 





= 1.956^^ 


+22 
- 4 


GeV 




= 1.9685t^GeV 


Md. 


= 2.002^^ 


+16 
-41 


GeV 


Md.o 


= 2.0067li GeV 




= 2.0821^ 


+22 
-31 


GeV 




= 2.1124t^GeV 



(47) 



For the lattice results, the first error is statistical and the second systematic. 
In this calculation the light quark in each meson is a normal quark. Thus the 
D* is the isospin-averaged vector state, with mass M^* = (Mj^.o + Md*±)/2. 
With somewhat large systematic errors, the largest of which comes from the 
lattice spacing, the spectrum is in broad agreement with experiment. 

It is also of interest to look at the My — Mp splitting as a function of the 
(pseudoscalar) meson mass. This can be extrapolated in inverse mass to the B 
scale, and to the infinite quark mass limit where the splitting should vanish. 
The splittings at both values of /? are given in Table 18 and are plotted for 
f] = 6.2 in Figure 14. There is little difference between linear and quadratic 
fits, so only the linear fit is shown. The splittings obtained from the linear fit at 
P = 6.2, using ro to set the scale, are taken as the central values quoted below. 
The second set of errors is systematic, obtained by combining in quadrature the 
differences arising from the following variations in procedure (in descending 
order of importance): using rUp rather than Tq to set the scale; (for the B 
meson) using a quadratic rather than linear extrapolation; multiple rather 
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Table 17 

The spectrum of heavy-light mesons with different definitions of the lattice spacing 
in GeV. 



state 




(3 = 6.2 


(3 = 6.0 




ro 


1.956^2 


1.952t^ 


Ds 


rup 


1.963^3 


1.958i^ 






1 A 

1.960+^ 


1 A 

1.958^^ 




ro 


2.002^6 


2.01lt^ 


D* 


nip 


1.9601:^ 


1.973l^ 




fn 


1.9731^15 


1.977tl^ 






2.082+^ 


2.0831^ 


D* 


nip 


2.052i^ 


2.053^3 






2.06X1:^ 


2.056lJ[! 




-50 ^ ' ' ' ' ' ' ' ^ 

0.2 0.4 0.6 0.8 



1/Mp(1/GeV) 

Fig. 14. The hyperfine splitting as a function of M in physical units at /? = 6.2. The 
circles have the scale set by ro and the diamonds by nip. 

than single exponential fits; using (3 = 6.0 data rather than f3 = 6.2 data. The 
results, with the comparison to experiment, are as follows: 

Lattice, this work Experiment [34] 

Md* - Md = ISQlllll MeV Md. - Md = 142.6 ± 0.5 MeV (48) 
Mb- - Mb = 2ltlt\l MeV Mb- - Mb = 45.2 ± 1.8 MeV 
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Table 18 

The hyperfine splitting of the heavy-light mesons at the B and D meson scale, in 
MeV. 





Uj 




— fi 




' u 


130+1 


140 +fi 


Mb* - Mb 


nip 




108 tl 




U 




111+14 






2ltl 


32 tl 


Mb* - Mb 


rup 


121? 


22 tl 




U 


1511 





6 Conclusions 

The decay constants of heavy-light and light mesons have been determined 
in the quenched approximation at two values of the coupling. The action 
and currents have been fully NP 0{a) improved. Good scahng of the decay 
constants is found. 

Uncertainties in the improvement and rcnormalisation coefficients are seen to 
have a large effect on the decay constants. The ALPHA and Bhattacharya et 
al. determinations of the improvement coefficients produce different values of 
/k, especially at /9 = 6.0. However, the /k do converge as the lattice spac- 
ing decreases. Computed quantities may well differ at 0{a?) at fixed lattice 
spacing but should agree in the continuum limit. 

The heavy-hght decay constants are found to be around 10% higher than, but 
compatible with, the world average [9] of quenched lattice determinations. 
HQS relations are found to be well satisfied by the data at both values of the 
coupling. The systematic uncertainties in the calculation have been discussed 
and, with the exception of quenching, estimated. The issue of 0{a?rn?) effects 
in the heavy quark extrapolations has been considered in various ways. In 
particular, the difference of the decay constants between (3 — 6.0 and /? = 6.2 
at fixed heavy quark mass, as a measure of discretisation effects, is small and 
remains so during the heavy quark extrapolation. The static point value of the 
decay constant, despite the lack of improvement, suggests the extrapolation is 
under control. Furthermore, HQS relations are satisfied and there is agreement 
between /bs//b computed by extrapolating the ratio and by taking the ratio 
of the extrapolated decay constants. 

The quenched spectrum of D mesons has been determined, shows good scaling, 
and is found to be in broad agreement with experimental data. 
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